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A structural model of hydrogenated amorphous silicon containing an isolated dangling bond is 
used to investigate the effects of electron interactions on the electronic level splittings, localization 
of charge and spin, and fluctuations in charge and spin. These properties are calculated with a 
recently developed density-matrix correlation-energy functional applied to a generalized Anderson 
Hamiltonian, consisting of tight-binding one-electron terms parametrizing hydrogenated amorphous 
silicon plus a local interaction term. The energy level splittings approach an asymptotic value for 
large values of the electron-interaction parameter U, and for physically relevant values of U are in 
the range 0.3 — 0.5 eV. The electron spin is highly localized on the central orbital of the dangling 
bond while the charge is spread over a larger region surrounding the dangling bond site. These 
results are consistent with known experimental data and previous density-functional calculations. 
The spin fluctuations are quite different from those obtained with unrestricted Hartree-Fock theory. 
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I. INTRODUCTION 

Amorphous silicon (a-Si) inevitably contains dangling 
bonds which lead to electronically active defect states 
in the band gap. For undoped material, most dangling- 
bond states are singly occupied, and their spins provide 
a well defined experimental signature. The Fermi level 
is controlled by the energy of the gap states. Hydro- 
genation of a-Si reduces the density of defect gap states 
by passivating the dangling bonds and thus restores the 
band gap, making hydrogenated amorphous silicon (a- 
Si:H) applicable to solar cell devices!!! However, even a 
small density of gap states can degrade performance, and 
gap states are also connected to degradation of device 
performance over time. Thus understanding the origin 
and properties of these states remains an important the- 
oretical challenge. 

The earliest theoretical work on defect states. Hf-iArSji 
and a-Si:H was based on tight-binding meAodsflMJolJ 
Biswas et al.0 and Fedders and CarlssonQ investigated 
the electronic structure of dangling and floating bonds in 
a-Si within tight-binding theory. They showed that the 
wave function of the gap defect states associated with 
the dangling bond is strongly localized on the threefold 
coordinated atomtl and relatively independent of strains, 
in contrast to the floating bond defect states. This differ- 
ence was taken to imply that the electron-spin resonance 
(ESR) signal in a-Si:H arises from dangling bonds. In 
tight-binding calculations without electron-electron in- 
teraction terms, the localization of the spin of the gap 
states is the same as that of the charge density, since 
the remaining occupied states do not adjust to the elec- 
tron charge in the dangling-bond gap state. This leads, 
in general to an overestimate of the charge density as- 
sociated with the gap state. More recently, density- 
functional calculations of dangling-bond states using the 
local-density approximation have been performed. a They 
yield a charge localization of less than 15% on the cen- 



tral atom. This finding at first appeared to be at variance 
with ESR experiments, which showed that over 50% of 
the spin density of the energy gap stajte-.is located on 
the central atom of the dangling bond.tULHl However, re- 
cent calculations using the local spin-density approxima- 
tion have shown that the degree of localization of the 
spin depsity is quite different from that of the charge 
densityJiil The energy cost to localize the charge den- 
sity is substantially larger than the energy to localize 
the spin densityMLj This demonstrates the importance 
of correlation effects for a correct description of the elec- 
tronic structure of the dangling bond, since in purely 
one-electron descriptions the charge and spin densities 
for a defect orbital are equivalent. 

Because of these correlation effects, the extent of 
the applicability of current implementations of density- 
functional theory to the electronic properties of defects 
in a-Si:H is not clear. These implementations break 
down in the limit of strong correlations. In addition, 
current density-functional codes do not provide informa- 
tion on the spin and charge fluctuations at the defect. 
For this reason, it is useful to study the defect states 
with a method that is valid in the limit of strong in- 
teractions, and that provides information on electronic 
fluctuations. In this work a recently developed method 
based on density-matrix functional theory, developed for 
isolated strongly interacting orbitalsliHl, is applied to the 
problem of a single dangling bond in a-Si:H. So far this 
method has only been applied to idealized models. In this 
paper we demonstrate that the method can be used to 
calculate the electronic structure of a semiquantitatively 
accurate model such as that treated here. Our results for 
the charge and spin distributions of the defect states are 
consistent with earlier results. We also make predictions 
for the fluctuations of the spin and charge. 

The paper is structured as follows. Section II describes 
the atomic structure of the model for a-Si:H. Section III 
introduces the Hamiltonian and describes the density- 
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matrix functional used to calculate the ground state en- 
ergy of the system. Section IV is the core of the pa- 
per, presenting our results for the electronic structure of 
the dangling bond. The effects of electron correlations 
on the energy of the defect state in the band gap are 
determined by comparison of results from the density- 
matrix fun ction al and the Hartree-Fock approximation 
in Section [VB. Charge and spin localization and fluc- 
tuations arc discussed in Section IV C. The results arc 



compared to density-functional calculations and experi- 



ments in Section [VD 



II. STRUCTURE MODEL OF 
HYDROGENATED AMORPHOUS SILICON 



The atomic model for amorphous hydrogenated silicon 
used here was employed in an earlier density-functional 
theory calculation.!!^ It contains 122 Si atoms and 20 hy- 
drogen atoms per fee unit cell, with periodic boundary 
conditions. The atomic positions were obtained by dou- 
bling a previous, smaller unit cell, and subsequently an- 
nealing the structure. The edge length of the fee cell is 
11 A. The hydrogen concentration of 14% is somewhat 
higher than what is commonly used in experimental sam- 
ples (c «10%). All hydrogen atoms are attached to the 
dangling bonds present in the structure, and each dan- 
gling bond is terminated by a hydrogen atom. 

To create a single dangling bond in the model of amor- 
phous hydrogenated silicon, hydrogen atom number 142 
was removed from silicon atom 108 and the structure 
was relaxed using a density-functional approach.E3 Fig- 
ure [j] shows the atomic structure surrounding the dan- 
gling bond site. It is seen that the dangling bond is 
surrounded by silicon atoms only, with hydrogen atoms 
further away. The orientation of the dangling bond is 
roughly in the [111] direction. 




III. HAMILTONIAN AND METHOD OF 
SOLUTION 



The electronic structure of the system is treated with 
a discrete Anderson-type Hamiltonian: 



Un ^n 0ll 



(1) 



where i and j are orbital indices, a denotes the spin and 
the dangling bond is orbital 0; cj^, Cjc, and riQ a are the 
usual second-quantized creation, annihilation, and num- 
ber operators. In the one-electron part of 7i the matrix 
elements hij include on-site energy terms and interatomic 
hopping terms. We restrict the interaction terms to the 
dangling-bond orbital because this orbital will have the 
greatest fluctuations in occupancy and will therefore be 
the most affected by the interactions. 

The angular dependences of the one-elcctroii|-terms 
are given by the Slater-Koster parameterization. til The 
Slater-Koster parameters are scaled with the interatomic 
distance, d, as 1/d 2 . There are several .tight-binding 
parametrizations in the literature for Si-H.aQ Knief and 
Niessen compared different tight-binding parameter sets 
for a-Si:H to the experimental density of states!] For a- 
Si:H they found a better agreement with experimental re- 
sults for the parameter set from Allan and MeleEl than for 
the tight-binding parameterization by Min et al.0 Both 
paramcterizations use orthogonal basis functions with a 
minimal basis set of s and p valence orbitals, and include 
nearest neighbor interactions only. In the following, the 
parameterization by Allan and Mele is used. 

The interactions of the up and down spin electrons of 
the dangling-bond orbital, described by the second term 
in 7i, were treated using a-recently developed density- 
matrix functional methodJla This approach treats cor- 
relations by including multiconfiguration effects in an 
approximate fashion. In previous tests for model ..sys- 
tems involving a single pair of interacting orbitalsE£l it 
was shown to give accurate results for several electronic 
properties for weak, intermediate, and strong electron- 
electron interactions. In the density-matrix functional 
method the ground state energy, (H), is approximated 
by a functional of the one-body density matrix, p, de- 
fined by pijc = {c\ a Cj a ). The expectation value of the 
one-electron part of the Hamiltonian is given exactly as 
a simple functional of the density matrix. The expecta- 
tion value of the interaction energy, 



FIG. 1: Structure of a-Si:H around the dangling bond site. In 
the left panel, the hydrogen in the center bonds to a Si atom, 
whereas in the right panel the hydrogen atom is removed pro- 
ducing a dangling bond. For clarity, the relaxations of the 
structure due to the removal of the hydrogen atom are not 
shown. 



Ei n t = U(n 01 n ol ), (2) 

is rigorously given as a functional of the local moments of 
the one-body density matrix projected on site O.Ea The 
exact form of this functional is not known. However, for 
systems with only two interacting orbitals, such as that 
studied here, a lower bound for the interaction energy 
holds which is given in terms of the second moment of 
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the density matrix: 

(Un 0T - £'int)(^(l - «ot - ™oi) + S int ) 
+ yjE int {Un ol -E int ). (3) 

A parallel result is obtained by switching up and down 
spins in the above inequality. In the "second-moment ap- 
proximation" that we employ here, the interaction energy 
is obtained as a function of p by replacing the inequality 
by an equality if this gives a positive value for the inter- 
action energy; if not the interaction energy is taken to be 
zero. 

The density matrix of the model system is taken to 
be that which minimizes the total energy, subject to 
the constraint that all of its eigenvalues must be be- 
tween zero and unity. This approach gives the correct 
density matrix for an exact density-matrix functionalt£l 
and is the appropriate avenue to use with our approx- 
imate functional. The resulting density matrix, unlike 
those obtained from density-functional calculations, has 
a range of eigenvalues between zero and one and is thus 
not idempotent (for U / 0). This is the correct behavior 
for interacting systems. The procedure for obtaining the 
energy-minimizing density matrix involves a constrained 
conjugate-gradient method described in more detail in 
Rcf. |l2j. The computer time required for the minimiza- 
tion is 0(N 3 ), where N is the number of orbitals. 



IV. ELECTRONIC STRUCTURE 
A. Without interactions 

Electronic densities of states of the model of the com- 
pletely hydrogenated amorphous silicon structure with 
U = were calculated using standard Brillouin-zone ip=. 
tegration techniques using a 5 x 5 x 5 k-point mesh.E£l 
The hydrogen passivatcs the dangling bonds and gives 
a well-defined gap of about 1.2 cV. Compared to exper- 
imental values of the energy gap pf 1.4 eV to 2.0 eV, 
the tight binding parameterizationB that we use under- 
estimates the gap. However, the density of states of our 
model of amorphous hydrogenated silicon compares well 
to the density of states for larger models such as the ones 
investigated by Holender and_Morgan using the same 
tight-binding paramcterizationO. 

In the single dangling bond electronic-structure calcu- 
lations, the one-body part of the Hamiltonian was trans- 
formed via the recursion methocO to a chain Hamilto- 
nian of length 80. Because the recursion steps quickly left 
the original unit cell, the cell was replicated periodically. 
The starting point for the recursion procedure was an sp 3 
hybrid orbital with the orientation of the above dangling 



bond orbital (see Sec. II). The chain was truncated at 
level 80 with no terminator. 

As a check on the accuracy of the recursion proce- 
dure, we obtain the electronic structure of a-Si:H with 
the dangling bond using both diagonalization of the 
chain Hamiltonian and standard Brillouin-zone integra- 
tion techniques. Figure || shows the density of states from 
the Brillouin-zone integration. The recursion method ob- 
tains a band gap of 1.4 eV, in comparison to 1.3 eV for 
the BZ method. The satisfactory agreement between 
these results indicates that the recursion chain is suffi- 
ciently long for an accurate description of the electronic 
structure. (The band gaps are different from those in 
the absence of the dangling bond because of the atomic 
relaxations and the finite size of the supercell). The 
dangling-bond gap state is found to be 73% localized on 
the dangling-bond orbital. The energy is 0.4 eV above 
the valence band edge and 0.9 eV below the conduction 
band edge. 




Energy [eV] 



FIG. 2: Electronic density of states of a-Si:H containing an 
isolated dangling bond calculated by Brilloin-zone integra- 
tion. The dangling bond leads to a defect state in the band 
gap. 



B. Effects of interactions on gap state energies 

To calculate the effects of interactions on the dangling- 
bond states, we apply the second-moment density-matrix 
algorithm described above to our Hamiltonian trans- 
formed into the chain representation. As U increases, 
the on-site energy of the dangling-bond orbital is low- 
ered by —U/2, keeping the average energy of the two 
interaction-split dangling bond states approximately con- 
stant. For comparison, we include results obtained by the 
unrestricted Hartree-Fock (UHF) method. 

The second-moment method as described above gives 
the total energy for a fixed number of electrons, but not 
directly the defect energy levels. In an interacting elec- 
tron system, a defect energy level is defined as a value of 
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FIG. 3: Energy levels of dangling-bond defect states in the 
band gap as a function of the interaction energy U. The en- 
ergy levels are given relative to the valence band edge obtained 
with the BZ integration method for [7 = 0. 
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FIG. 4: Spin of the dangling-bond orbital as a function of 
interaction energy U for the case that the chemical potential 
lies between the two defect levels. The spin is measured in 
units of the electron spin. 



the chemical potential \x at which the number of electrons 
in the system changes abruptly. Minimizing the thermo- 
dynamic potential E(N)—fxN (at zero temperature), one 
readily shows that the number of electrons changes from 
N to N + 1 when 



H = E(N + 1)-E(N), (4) 

which is taken to be the gap state energy. The valence 
and conduction-band edges are defined in a similar fash- 
ion (they are discrete states because our chain has fi- 
nite length) . Figure || shows the gap energy levels and 
band edges obtained in this fashion as functions of the 
interaction energy, U, for the second-moment and the 
UHF approximations. We first note that the conduc- 
tion and valence-band edges depend only weakly on the 
Coulomb repulsion on the dangling bond. The energy 
of the gap states in the second-moment approximation 
varies roughly linearly with U, for small U, as expected 
on the basis of first-order perturbation theory. At larger 
values of U, beyond about 3 eV, the splitting approaches 
a finite limit. The UHF results are in agreement with 
the second-moment results up to about U = 2 eV, but 
the levels continue to split linearly with energy until they 
merge with the valence and conduction bands. On the 
basis of exact diagonalization many-body calculations for 
small clusters including only the nearest few orbitalsEJ, 
we feci that the behavior of the second-moment approx- 
imation is correct. In the diagonalization calculations 
we find that it is possible to add a second electron to 
the gap states without an energy increase proportional 
to U, because the inclusion of correlation effects allows 
the electrons to avoid both being in orbital at the same 
time. 



C. Spin and charge on dangling-bond site 

The permanent spin on the dangling-bond site is ob- 
tained directly in terms of the appropriate elements of 
the density matrix: 

(s ) = (not -noi). (5) 

Figure |I] shows the dependence of this spin on U for the 
case when the chemical potential is between the two de- 
fect levels. In the second-moment approximation, the 
spin rises for small values of U and then levels off at a 
value about 15% higher than the zero-t/ value. In con- 
trast, the spin in the UHF approximation continues to 
rise at the highest values of U that were treated, and 
eventually approaches unity. This is analogous to the be- 
havior observed for Anderson-chain modelsE3. In these 
models, the UHF approximation overestimates the lo- 
cal moment on the dangling-bond site in order to reduce 
the interaction energy. The second-moment approxima- 
tion, however, does not yield such a large local moment. 
The reason is that in the second-moment approximation 
correlations are included via multiconfiguration effects, 
rather than by varying the moment of a single configu- 
ration. 

The fluctuations in the spin and charge are also ob- 
tained straightforwardly by the density-matrix approach. 
To accomplish this, we note that Eq. @ implies that 
E-mt/U = {no^riQi). Then the spin and charge fluctua- 
tions on site are obtained as 

<(A So ) 2 } = (s 2 ) (so) 2 (6) 

= (no? + "o;> - ("or - "oi) 2 - 2E int /U 

and 

((An ) 2 ) = {nl) (n ) 2 (7) 

= (n | + n Q \) - (not + ™o|.) 2 + 2E iat /U. 
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FIG. 5: Fluctuations of the spin (a) and the charge (b) on the 
dangling-bond orbital as a function of the interaction energy 
U. 
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FIG. 6: Spin and charge density of the defect state projected 
on the chain sites for U — 3 eV. 



Here we have used the fact that (HqA = (not) an d 
( n oi) = ( n oi)- The dependence of ((As ) 2 ) on U is plot- 
ted in Figure |^(a). In the second-moment approximation, 
the spin fluctuations increase with U , to an asymptotic 
value about 60% higher than the U — value. In con- 
trast, the UHF results reveal a monotonic decrease of 
((Asfj) 2 } with U. These-.results are consistent with the 
Anderson-chain rcsultsJ13 The decrease in the spin fluc- 
tuation with U in the UHF approximation results from 
the increased moment obtained in this approximation, 
while the increase observed in the second-moment ap- 
proximation results from the reduced occupancy of the 
zero-spin states in which both orbitals on site are empty 
or filled. The charge fluctuations on site are shown 
in Figure |E|(b). In both the second- moment and UHF 
approximations, the fluctuations drop with increasing 
U. However, the second-moment approximation yields 
a more pronounced drop than the UHF approximation. 
The behavior of the charge fluctuations, like that of the 
spin fluctuations, is due to the suppression of configura- 
tions with zero or double occupancy with increasing U. 



D. Comparison to Experiments and Previous 
Theory 

The primary outputs of these calculations that can be 
compared with experiment and previous theory are the 
spin and charge localization of the gap states, and their 
splitting. The components of the gap-state charge den- 



sity on the sites other than arc found by evaluating 
the changes in site-projected charges when the chemical 
potential crosses up through the lower gap level. The pro- 
jected charge density on the sites of the recursion chain 
is shown in Figure H(a) for the second-moment approxi- 
mation and the UHF approximation for an intermediate 
value of the Coulomb repulsion, U — 3 eV. The neglect 
of the correlation effects in the UHF approximation re- 
sults in an overestimate of the charge localization of the 
defect state on the dangling bond orbital. In the second- 
moment approximation, the charge is strongly delocal- 
ized over a large part of the chain. The spin on the other 
hand is strongly localized on the dangling bond orbital. 
The UHF yields a slightly larger spin than the second- 
moment approximation. Overall, the correlation effects 
lead to a larger degree of localization for the spin than for 
the charge. This confirms the LSDA results by Fedders 
et alfi 

We expect the electronic structure of gap states in a- 
Si:H to be described qualitatively by the present model. 
A quantitative comparison to experimental data is diffi- 
cult since the strength of the bare Coulomb interaction, 
U, for dangling bonds in amorphous silicon is not known 
precisely, and we ignore relaxation processes, which are 
known to reduce the effective correlation energy.Ej How^. 
ever, most estimates of U are in the range 2-5 eVESEHl 
In this range, the gap state splitting is already close to 
its asymptotic large-?/ value, which is determined by the 
coupling of the gap state to the neighboring orbitals. For 
this reason, rather than presenting a single set of results, 
we present a range of results corresponding to the above 
range of values of U. These are compared to experimental 
data and previous density-functional results in Table Q. 
For completeness the UHF results for a smaller range of 
Coulomb interactions, C/ = 2 ... 3 eV, are also included 
in the comparison. Beyond these values of U the defect 
states are no longer in the gap in the UHF approxima- 
tion. 



TABLE I: Comparison of the splitting Ae of the defect state in 
the energy gap as well as the spin and charge localization from 
different methods to experimental values. The projected spin 
(sdb) and charge (nab) of the defect onto the four sp 3 orbitals 
of the atom associated with the dangling bond are given in 
units of the electron spin and charge respectively. 





Ae [eV] 


(Sdb) 


(n d b) 


LDA&-, 






0.10-0.15 


LSDAB 


0.25-0.30 


0.41-0.52 


0.16 


UHF 


0.6-0.9 


0.65-0.70 


0.40-0.44 


2 nd -moment approx. 




0.60-0.S3-, 


0.10-0.29 


Experiment 


oi-o.fil 


0.50-0.80eIH 











ESR and photolumincsccncc spectroscopy measure- 
ments have given values for the splittrng-pf the two gap 
states ranging from 0.3 eV to 0.4 eVtiJo. These split- 
tings are close to those obtained here for a wide range of 
values of U. Comparable agreement is obtained by the 
LSDA calculations, but the UHF method substantially 
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overestimates the splittings. The extent of spin localiza- 
tion in the present results is very insensitive to U, and is 
roughliiiiJ-thc middle of the range obtained in ESR exper- 
imentsBBIill Again, the LSDA results are quite compa- 
rable. In both approaches, the degree of spin localization 
of the defect state on the dangling bond is much greater 
than the localization of the charge; however, this does 
not hold for the UHF results. Overall, the agreement 
of the results of the second-moment approximation with 
experimental values is surprisingly good, considering the 
simplicity of the underlying tight-binding model. To our 
knowledge, no experimental methods exist for measuring 
the extent of charge localization on the dangling-bond 
orbital. 



V. CONCLUSION 

The above results illustrate the applicability of the 
second-moment implementation of density-matrix func- 
tional theory to electronic-structure models with semi- 
quantitative accuracy such as the tight-binding model 
used here. The results show that the splitting of the 
gap states is smaller than expected from Hartree-Fock 
calculations and approaches a finite limit for large val- 
ues of the Coulomb repulsion. This effect can by ex- 
plained by the enhanced correlation of the electrons in 
dangling-bond states with increased Coulomb repulsion. 



It is found that the spin of the defect state is strongly 
localized on the dangling bond orbital while the charge 
is quite delocalized. These results are rather insensitive 
to the specific value of the Coulomb repulsion parame- 
ter, and are in fairly good agreement with results from 
electron spin resonance experiments and local-spin den- 
sity functional calculations. Our results for the charge 
fluctuations are similar to those obtained from Hartrec- 
Fock theory, while the results for spin fluctuations are 
quite distinct. We are not aware of existing methodolo- 
gies for measuring these fluctuations, but such measure- 
ments could provide an accurate test of the precision of 
the methods used here. ,— ■ 

Because of the previously demonstrated^ applicabil- 
ity of the second-moment implementation of density- 
matrix functional theory to strongly interacting systems, 
it would be desirable to apply it to transition-metal im- 
purities in both semiconductors and insulators. At this 
point, such applications cannot be performed because we 
do not have a suitable energy functional for such a multi- 
orbital impurity. Future work in this field should aim to 
extend the present methodology to include such systems 
with more than two interacting orbitals. 
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